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CN . Abstract 



(^T) • We explore interpretations of the anomaly observed by HI and ZEUS at HERA in 

deep- inelastic e + p scattering at very large Q 2 , in terms of possible physics beyond the 
Standard Model. Since the present data could be compatible with either a continuum 



{NJ ■ or a resonant solution, we discuss both the possibilities of new effective interactions and 

the production of a narrow state of mass M ~ 200 GeV with leptoquark couplings. 
We compare these models with the measured Q 2 distributions: for the contact terms, 
Q^ ■ constraints from LEP 2 and the Tevatron allow only a few choices of helicity and flavour 

structure that could roughly fit the HERA data. The data are instead quite consistent 
Q-i! with the Q 2 distribution expected from a leptoquark state. We study the production 

cross sections of such a particle at the Tevatron and at HERA, the latter in the cases 
where it is produced from either a valence or a sea quark. The absence of a signal at the 
Tevatron disfavours the likelihood that any such leptoquark decays only into e + q. We 
then focus on the possibility that the leptoquark is a squark with i?-violating couplings. 
^ In view of the present experimental limits on such couplings, the most likely production 

channels are e + d — > cl or perhaps e + d —* t, with e + s — > t a more marginal possibility. 
We point out that the cl could have competing branching ratios for i?-conserving and 
i?-violating decay channels, whereas t decays would be more likely to be dominated by one 
or the other. Possible tests of our preferred model include the absence both of analogous 
events in e~p collisions and of charged current events, and the presence of detectable 
cascade decays whose kinematical signatures we discuss. This model could also make an 
observable contribution to K — > txvv and/or neutrinoless (3/3 decay. We also discuss the 
possible implications for the Tevatron and for e + e~ — * qq and neutralinos at LEP 2. 
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1 Introduction and Summary 



The HERA experiments HI |TJ and ZEUS have recently reported an excess of deep-inelastic 
e + p scattering events at large values of Q 2 > 1.5 x 10 4 GeV 2 , in a domain not previously 
explored by other experiments. With a total e + p integrated luminosity of 14 pb _1 , HI [l|] 
observes 7 events with large e + -jet invariant masses M = y/xs, clustered around M = 200 
GeV, in which the positron is backscattered at large y = Q 2 /M 2 . Similarly, ZEUS || with an 
integrated luminosity of 20 pb _1 observes 5 events at comparable large values of Q 2 , x and y. 
Although the HI and ZEUS data are mutually consistent and the presence of the same type 
of excess in the two experiments is certainly impressive, the detailed features of the events are 
not exactly the same in HI and ZEUS. The events of HI are more suggestive of a resonance 
with e + -quark quantum numbers than the ZEUS data points, which are more scattered in 
mass. The difference could, however, be due to the different methods of mass reconstruction 
used by the two experiments, or to fluctuations in the event characteristics. Of course, at this 
stage, due to the limited statistics, one cannot exclude the possibility that the whole effect is a 
statistical fluctuation. This will hopefully be clarified soon by the coming 1997 run. Meanwhile, 
it is important to explore possible interpretations of the signal, in particular with the aim of 
identifying additional signatures that might eventually be able to discriminate between different 
explanations of the reported excess. 

Since the observed excess is with respect to the Standard Model expectation based on 
the QCD-improved parton model, the first question is whether the effect could be explained 
by some inadequacy of the conventional analysis without invoking new physics beyond the 
Standard Model. In the case of the apparent excess of jet production at large transverse energy 
Et recently observed by the CDF collaboration at the Tevatron ||, it has been argued |4|] 
that a substantial decrease in the discrepancy can be obtained by modifying the gluon parton 
density at large values of x where it has not been measured directly. In the HERA case |], 3 , 
a similar explanation is apparently not viable. In this case, the valence quark densities are the 
most relevant ones, and they have been measured directly |j] in the same range of x at much 
lower values of Q 2 . Since the values of x which are relevant for the HERA excess are quite 
large (x ~ 0.5), it is possible in principle that higher-order effects of the Sudakov type, not 
accounted for by the standard next-to-leading-order QCD analysis of the structure function 
data, could affect the low-energy extraction of the partonic densities and their evolution to 
high Q 2 ||. It should be remarked, however, that the most recent measurements of F2 from 
the HERA experiments explore the same high-x values in a large range of Q 2 values, and no 
anomaly was found at large x up to Q 2 ~ 10 4 GeV 2 [p], [|. The evolution logarithms could not 
explain the abrupt occurrence of the effect, which is undetected at Q 2 ~ 10 4 GeV 2 but fully 
visible at Q 2 ~ 2 x 10 4 GeV 2 . Therefore one can safely conclude that Sudakov effects cannot 
provide a credible explanation of the observed excess. We also note that, if the parton densities 
were to blame, very similar effects should be seen in both neutral and charged current channels, 
with both e + and e~ beams. This can be checked in the near future. We do not consider this 
alternative in the following, but concentrate on interpretations based on possible new physics. 

We first discuss the possibility that the observed excess is a non-resonant continuum. Within 
this scenario, a rather general approach is to interpret the HERA excess as due to an effective 
four-fermion eeqq contact interaction [|7| with a scale A of order 1.5-2.5 TeV. It is interesting 
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that a similar contact term of the qqqq type, with a scale of exactly the same order of mag- 
nitude, could also reproduce the CDF excess in jet production at large E T We study the 
contact interaction scenario for the HERA excess in some detail. In order to interfere with the 
electroweak gauge interactions, the contact term is taken as the local product of two mixed 
vector- and axial-vector currents. We study the x and Q 2 distributions that correspond to 
different flavours, signs and choices of helicity for the contact terms (i.e., LL, RR, LR or RL), 
and compare them with the HERA data. Strong bounds on the possible magnitudes of the 
interaction scales A are imposed by CDF data || on the Drell-Yan production of e + e~ pairs 
at large invariant mass, and by LEP 2 data on hadron production |L0], [TTj] . If we restrict our 
analysis to one particular term at a time, though in general one cannot exclude a superposition 
of different chiral structures, we find that most of the individual contact terms that could fit 
the HERA data are already excluded. Only for particular choices of quark flavour, sign and 
helicity can one obtain even rough agreement with the HERA data while escaping the existing 
bounds. We present examples of these models, and point out the desirability of further tests 
at LEP 2, where a complete analysis by all the experimental collaborations is still lacking, and 
at the Tevatron. It is interesting that the existence of the appropriate contact terms could be 
soon excluded, or their effects discovered in these experiments. We recall that the effects of 
contact terms should be present in both the e + and the e~ cases with the same intensity, and 
possibly also in the charged current channel, if left-handed currents are invoked. 

We then focus on the possibility of a resonance with e + q quantum numbers, namely a 
leptoquark. Most probably the production at HERA occurs from valence u or d quarks, since 
otherwise the coupling would need to be quite large, and more difficult to reconcile with existing 



limits (12], [13], [XT], JJj. Assuming an S'-wave state, one may have either a scalar or a vector 
leptoquark. Although we mostly consider the first option, we also include some discussion of 
the vector case. Defining the coupling A for a scalar (f> by X^e^q^ or A^e^g^, the observed excess 
of ~ 10 events in 34 pb _1 , observed with an efficiency of ~ 80%, suggests values of A ~ 0.025 
or 0.04 for production from u or d quarks respectively. The corresponding natural decay 
width is of the order of a few MeV. We compute the Q 2 distribution predicted by leptoquark 
production, and show that it matches the data better than the corresponding distributions for 
contact terms. A scalar e + u or e + d state couples to the following SU(2) doublet combinations: 
e^iuL.di), (e^, z/r)wr, or (e^^RjdR. This implies that a scalar with e + u or e + d quantum 
numbers can decay into vq final states only if it has another Yukawa interaction besides that 
responsible for its production. This additional interaction involves a lepton field of opposite 
chirality and is strongly constrained by pion decays. In its absence, a leptoquark would not 
be able to explain any resonant signal in the charged- current channel. The HI collaboration 
has reported |IJ four events in this channel, with a Standard Model background of about two, 
but ZEUS @ has not reported a recent charged- current analysis. We note that the situation is 
different for a vector leptoquark in the e + q case, or for a scalar leptoquark in the e~q case. In 
the latter case we could, for example, have a coupling to the weak isospin singlet e~[UL — Vhdh 
that indeed leads to both neutral and charged current decay modes. 



Leptoquarks would be produced via QCD interactions at the Tevatron [|15|]. We find that 
a scalar leptoquark of mass M ~ 200 GeV has a production cross section of around 0.2 pb 



4 Note, however, that this interpretation is not strengthened by more recent data on the dijet angular distri- 
bution 0. 
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at the Tevatron. The cross section for vector leptoquarks is somewhat model-dependent, but 
expected to be much larger, as discussed later. Given the large value of the cross section, a 
leptoquark branching ratio B(e + q) < 1 into the observed e-jet channel is perhaps needed, even 
in the scalar case, to avoid a possible combined CDF/DO exclusion limit. The current best 
limit from the Tevatron for a first-generation leptoquark is 194 GeV for B(e + q) = 1, recently 
given by the DO collaboration [ i"6fl , and the corresponding limit for B(e + q) = B{uq) = 0.5 is 



143 GeV [16| . Thus any scalar with leptoquark quantum numbers might need additional decay 
modes beyond those given by the A interaction introduced above for its production mechanism. 
Perhaps the most appealing form of leptoquark is a squark Jl7| with couplings that violate 



R parity |18fl . This possibility has been put forward in connection with the HERA events also 



in ref. |19[ . In terms of supersymmetric chiral multiplets, the relevant coupling is given by 

KikUQjDl (l) 

where L iy Qj and D% are superfields of lepton doublets, quark doublets and quark singlets 
respectively, and i,j, k are generation indices. Leaving production from the sea aside for the 
moment, the processes relevant for HERA that arise from this coupling are e^dR — > ul or cl 
or ti- We find that the first possibility is eliminated by existing limits on A' m , in particular 
those from (3(3 decay f2~0fl , whereas the latter are still permitted. In the following, we study the 



scharm c and stop t possibilities in some detail. We recall that the i?-violating decays of the 
produced squark must compete with the ordinary i?-conserving decays. Although some sizeable 
additional decay channels are welcome in view of the non-observation of a signal at CDF/DO, 
the i?-conserving channels should have a moderate rate, otherwise the coupling A- jfc required 
to explain the HERA excess becomes dangerously large, particularly in view of upper limits on 



A' 121 coming from (3(3 [20| and K — > nvv (gTJ decays. We make a careful study of the regions 
of parameter space for the supersymmetric model where the balance of /^-conserving and R- 
violating decays is favourable, for both the c and t cases. Such a balance is more likely in the 
cl case than for the t. In the squark scenario there would be no signal in the charged-current 
vq channel, but there could be .R-conserving decay signatures, whose kinematical properties 
we discuss later. No signal is expected from the same sparticle with e~ beams, unless one is 
sensitive to production from the d sea density. A distinctive signature of the c possibility could 
be the appearance of a signal in K — > -nvv close to the present upper limit. 

We have also examined possible signatures of i?-violating supersymmetric interactions at 
LEP 2. We find that interference effects in the reactions e + e~ — > qq are unlikely to be detectable, 
unless the HERA squark is produced from the sea, in which case a signal might be detectable 
in ss final states. However, other effects of i?-parity violation could be observable, such as 
e + e~ — > where x° denotes the lightest neutralino, thanks to its i?-violating decays. We 

also discuss the compatibility of the squark explanation of the HERA events with the model 



recently proposed |22| to interpret the four-jet anomaly found by ALEPH ||23|| , but not seen by 
the other LEP experiments 
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Table 1: Preliminary OPAL 95% CL limits on the effective contact interaction 
scale Afj (in TeV). 



2 Effective Contact Interactions 

Whereas the HI data are at first sight quite suggestive of the production of a resonance in 
the s channel, the spread in x 2 a of the ZEUS data seem to favour the possibility of an 
effective 4-fermion contact interaction, which we pursue first as a more conservative option. As 
is customary in the literature we parametrize the contact interactions in terms of the mass 
scale A appearing in the following effective Lagrangian: 

A = 4tt -UTvl e^ e i Qjl^j • ( 2 ) 

i,j=L,R VHj) 
q=u,d 

We allow for independent couplings of u and d quarks, as well as for independent couplings of 
all different helicity states Q. The parameter rjij takes the values ±1, and allows for constructive 
and destructive interferences in the different channels. 

Very tight constraints on the size of such possible interactions have been set in the past |25|. 



Recent preliminary results from dielectron production at the Tevatron || and from hadron 
production at LEP 2 ||T0| , |TT| restrict even further the allowed ranges of the parameters A?- in 
(0). The most recent analysis by OPAL |11|| , in particular, sets the 95% CL limits on the 16 
independent couplings in @ shown in Table |l| The CDF limits || have only been given for 
the isoscalar u + d quark flavour combination, and for the LL helicity combination. We have 
simulated the effects in hadronic collisions of the contact interactions for which the OPAL limits 
allow good fits of the HERA data. For these cases we extracted from the Tevatron dilepton 
data limits on the relevant parameters A^ by analogy with the limits provided by CDF for the 
isoscalar LL case 0. 

Since contact interactions do not generate any particular structure in the x distributions, 
we discuss their impact on the integrated Q 2 distributions of the HERA data, as provided in 
their papers and combined in their public presentations ]I], . We have calculated the effects of 

5 The variable xi a is the Bjorken x variable extracted using the double- angle method ||, |). We discuss in 
the Appendix issues related to this spread and to the comparison between values of x reported by HI H and 
ZEUS H , in the light of possible initial-state radiation. 

6 We do not consider here scalar current couplings, because they are very strongly constrained by low-energy 
data on helicity-suppressed decays . 
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Figure 1: HERA data for the integrated Q 2 distribution, compared to current limits 
on effective contact interactions from OPAL and CDF. Only the four combinations 
that best reproduce the HERA data are shown. The lower solid line corresponds to the 
prediction of the Standard Model. The upper solid curve corresponds to the decay of a 
200 GeV s- channel resonance, produced at a rate compatible with the reported excess 
of events. 



each one of the 16 4-fermion couplings, including its interference with the Standard Model DIS 
processes. In the cases where the OPAL constraint allows a fit, we have applied the inferred 
CDF constraints on the corresponding AfL-, which are usually stronger. 

We applied the analysis cuts of the HI and ZEUS experiments, as described in their publica- 
tions, and combined the expectations for the respective integrated luminosities and efficiencies []. 
Only a few of the 16 possible couplings are at all compatible with the HERA data, and the four 
best cases are presented in Fig. [1|. In none of these cases is the agreement in shape between 
data and expectations particularly good. We note the essential roles played by both OPAL and 
CDF in constraining the possible effective contact interactions: in particular, there are good 
fits of the HERA data for couplings with magnitudes that are at best compatible with the 
OPAL data alone, such as the choice K\ L = 1.4 TeV, with rj = — 1. However, this possibility 
is excluded by the CDF limit, which is stronger for this specific coupling, namely larger than 

7 We have verified that our Standard Model predictions, after accounting for the analysis cuts, efficiencies and 
integrated luminosities of the two experiments, agree with those presented in the HI and ZEUS papers [H, ||. 
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3 TeV. 

We conclude from this study that, while the contact interaction hypothesis cannot be entirely 
ruled out, the strong constraints already set by the LEP and Tevatron experiments do not allow 
for good fits of either the event rates or the Q 2 distributions of the HI and ZEUS data. We 
remark once more that, in any case, for the contact interaction hypothesis to be tenable, the 
apparent resonant structure in the invariant e + q mass distribution reported by HI could only 
be the result of a statistical fluctuation, and should be washed out by higher statistics. We 
would also expect that a joint effort of the four LEP collaborations using the combined set of 
all LEP 2 data already on tape could further restrict the allowed ranges of the A parameters, 
as could a combined analysis of CDF and DO data. 



3 Leptoquarks 

Since the excess found by HI [1|] occurs in a small range of e + q invariant masses 0, it is natural 
to examine models containing a new boson with leptoquark quantum numbers, which may 



be classified according to their spin and isospin quantum numbers ||26|1 . As discussed in the 
introduction, they may couple to a fermionic current constructed out of a quark and a lepton 
with a coupling constant A. In the narrow-width approximation, the leading-order parton-level 
cross section for production of a leptoquark of spin J at HERA is given by: 

a = f X 2 Fj (F = l, F x = 2). (3) 
The convolution of the parton cross section (|J) with the parton densities of the proton [27 



and with the effects of initial-state photon radiation from the positron, yields the cross sections 
shown in Fig. |2|. We have assumed in this figure a leptoquark coupling to only one fermionic he- 
licity, and have used A = 0.01 as a reference value. Figure |2| shows results for all quark flavours. 
As is clear from the figure, all effective sea quark luminosities are significantly suppressed for 
masses in the 200 GeV region, where the hypothetical HERA signal lies. The cross sections for 
A = 0.01 in the different production channels to produce a leptoquark with mass 200 GeV are 
given in Table [2[ We have verified that the acceptance cuts imposed by HI and ZEUS on 
their data have efficiencies of approximately 80% and 100%, respectively. Accounting for the 
detector efficiencies of the two experiments as quoted in their papers (of the order of 80%), for 
the relevant integrated luminosities, and assuming a total of 10 signal events in the combined 
experiments, we find that a value of A of approximately 0.04/v^B is required for leptoquark 
production by e + d collisions, with correspondingly larger couplings required for production by 
e + collisions with sea quarks. The cross sections for a vector leptoquark are a factor of two 
larger, and hence would require couplings smaller by a factor of V2. The implications of the 
value A ~ 0.04 in the case of scalar quark production via an i?-violating interaction will be 
discussed in the following section. 

8 Although the ZEUS events are more spread in invariant mass, and appear at larger values of X2 ai we note 
that an inter-collaboration working group has found that the two data sets are compatible, and that ISR and 
other instrumental effects could in principle cause shifts of several % in the observed values of x^a , as discussed 
in the Appendix. 
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Figure 2: Scalar leptoquark (or R-violating squark) production cross sections at 
HERA, including the effects of initial-state radiation (ISR). The contributions of dif- 
ferent quark flavours in the proton are shown separately. 



o-(fb), A = 0.01 


e + u e + d e + d e + u e + s e + c e + b 


(with ISR) 
(no ISR) 


106 25.6 0.98 0.47 0.47 0.23 0.12 
117 28.4 1.12 0.53 0.53 0.26 0.14 



Table 2: Production cross sections for a scalar leptoquark (or R-violating squark) of 
mass 200 GeV at HERA, assuming a coupling A = 0.01, showing the effects of ISR. 
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It is important to explore the possible implications of the existence of 200 GeV leptoquarks 
at the Tevatron, since there leptoquarks are produced via model-independent QCD processes 
with potentially large rates. The production cross section for scalar, colour-triplet particles 
at the Tevatron is given as a function of the mass in Fig. |3], using the MRSA' pTJ parton 
distribution function set, and a renormalization/factorization scale /i = m. The lowest-order 
production cross section for a pair of 200 GeV leptoquarks is 0.18 pb 0. The total integrated 
luminosity collected by the CDF and DO experiments is of the order of 200 pb -1 , which would 
yield approximately 36 events if the hypothetical leptoquarks had a 100% branching ratio into 
electrons. The detection efficiencies for such a signal quoted by the CDF and DO experiments 
are each of the order of 20%, resulting in about 7 detected events in total. Although the 
expected 3.5 events per experiment are not sufficient to exclude a leptoquark of 200 GeV in 
either CDF or DO individually, the best current limit being 194 GeV by DO as mentioned in the 
Introduction, it is likely that exclusion of a 200 GeV scalar leptoquark with 100% branching 
ratio into electrons at more than the 95% CL would result from the absence of a signal in a 
combined analysis of the data collected by CDF and DO. For this reason, in the following section 
we disfavour models with decay modes dominated by electrons although this possibility is 
not yet rigorously excluded. 

The case of a 200 GeV vector leptoquark is most likely totally ruled out by the Tevatron data, 
since the production rate can be as much as a factor of 10 larger than that of scalar leptoquarks, 
as shown in Fig. [| We emphasize that in calculating this curve we have only included the light 
quark annihilation processes, and assumed minimal couplings to the gluons [^J . In the absence 
of a definite model of vector leptoquarks, their coupling to gluons is not uniquely defined, and 
in general leads to bad high-energy behaviour and unphysically large cross sections. This 
problem has been studied in detail in refs. [^, |3T|, where J = 1 leptoquark-pair production 
was considered for a general class of anomalous couplings. In ref. [[30], |31| it was shown that, 
even allowing for anomalous couplings, and selecting them so as to minimize the production 
cross section via destructive interference, the total rate would still be a factor of two larger 
than that for scalar leptoquarks. In the absence of a signal, such a large rate would not be 
consistent with a combined CDF+D0 analysis if B(e + q) = 1. If one discards the possibility of 
such a fine-tuning in the anomalous couplings, values of B(e + q) significantly below 1 could be 
excluded. 



4 R- Parity Violation 

We find it attractive to embed a hypothetical leptoquark in a well-motivated theoretical frame- 
work capable of constraining its properties and providing other experimental signatures, as is 

9 For our choice of renormalization scale, we also expect a multiplicative K factor of the order of 1.10 to 1.15 
due to higher-order QCD corrections. These have been evaluated in the case of supersymmetric scalar quarks 
in psj . The case of leptoquarks can be recovered by assuming a very large gluino mass. In this case, diagrams 
due to four-squark operators which are not present in the leptoquark case only appear in the gluon-fusion 
production channel, which is significantly suppressed. 

io p r i or ^ generic leptoquarks could avoid this problem by having two different Yukawa couplings, allowing 
vq decays as mentioned in the Introduction, or by coupling to fj,q or rq. However, the fiq option is strongly 
constrained by limits on flavour-changing interactions, such as fi — e conversion on nuclei |l^, |i~3|| . 
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Figure 3: Scalar and vector leptoquark production cross sections at the Tevatron. In 
the case of the vector leptoquark, as discussed in the text, we have only included the 
light qq annihilation processes. 
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provided by the minimal supersymmetric extension of the Standard Model [17| with violation 
of R parity ||18|| . The corresponding superpotential can be written in the form 



W R = tuff Li + XiikULiEi + X'^.L.QjI)', + X'^DPt , (4) 

where H, L i: Ej,Qk, (U, D)i denote superfields for the Y = 1/2 Higgs doublet, left-handed 
lepton doublets, lepton singlets, left-handed quark doublets and quark singlets, respectively. 
The indices i,j, k label the three generations of quarks and leptons. Henceforth, we work 
in a basis for the Lj and the Y = —1/2 Higgs doublet H in which /z, = 0, and the only 
surviving bilinear term is the Higgs mixing fiHH. Furthermore, we assume the absence of the 
A" couplings, so as to avoid rapid baryon decay, and the A couplings play no role in our analysis. 



4.1 Production mechanisms 

The squark production mechanisms permitted by the A' couplings in ([|) include e + d collisions 
to form ul,cl or t i; which involve valence d quarks, and various collisions of the types e + di 
{i = 2, 3) or e + Ui {i = 1, 2, 3) which involve sea quarks. The required magnitude of the coupling 
A' is fixed by the observed product of the cross section a and the squark branching ratio B for 
the .R-parity violating mode q —>■ e + q'. From the results of the previous section, summarized 
in Table 2, we infer that the valence production mechanism requires A^ (j = 1,2,3) to be 
about 0.04/ y/B, while any of the sea production mechanisms require A' ljA , > 0.3/ \f~B (j,k = 
1,2,3). The latter are only marginally compatible with LEP 2 limits [11], and with previous 
HI limits ]14] in the cases j or k = 1. 

The required values of the A' ljfc are to be compared with the upper limits available from 
various other laboratory experiments 0. It has been inferred from upper limits on neutrinoless 
(3/3 decay that J20J 

where rrig is the mass of the lighter of ul and da, and rrig is the gluino mass. This limit 
excludes any production mechanism involving only first-generation particles, and in particular 
the valence parton process e + d — ► ul- 

For charm squark production e + d — > Cl, the most important constraint on the relevant 
coupling constant A' 121 comes from limits on flavour-changing neutral current processes. The 
simultaneous presence of several A' couplings with different flavour indices leads in general to 
dangerous tree-level flavour violations. Usually one makes the most conservative assumption 
that only a single A' coupling with specific flavour indices is non-negligible. However, because 
of the mismatch in flavour space between the up- and down-type left-handed quarks, this 
hypothesis cannot be simultaneously satisfied both in the up and down sectors. We will work 
in a basis where, in terms of the lepton (Ni,Ei) and quark (Ui,Di) mass eigenstates, the A' 
interaction term in the superpotential is written as 

A^(iV^A - E^Dl , (6) 



11 We will not consider here very stringent limits on i?-parity violating interactions coming from cosmological 
considerations of the baryogenesis energy scale ]3^| , since there are ways to avoid these in principle, such as 
baryogenesis at the electroweak scale |33l . 
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where Vij is the usual Cabibbo-Kobayashi-Maskawa matrix. We will also implicitly assume 
that the only sources of flavour violations are described by V and by the .R-parity violating 
interactions. Because of the non-trivial mixing in the down sector, the A^ couplings are 
bounded by B(K^ 



7T 1 VV) 



< 2.4 x 1(T 9 to be ET1 



\K jk \ < 2 x 10- 



m 



200 GeV . 



for j = 1,2, k = 1,2,3. 



(7) 



Therefore a cl interpretation of the HERA data, which implies A' 121 

400 GeV 

m d R > 



0.04/vB, is possible if 



However, this bound on vng R can be partially relaxed if the mixing in the down sector is 
somewhat suppressed by the simultaneous presence of various non-vanishing coupling constants 
X' ijk . For instance, allowing for several X[ x with different indices j, the bound in eq. (0) becomes 



\ 
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V-2 



21 



V A' Vn 

I V22 



< 2 x 10" 2 



m 



200 GeV 



(9) 



If A' m saturates the bound in eq. ((^), then A' 121 can be as large as 4 x 10~ 2 , even with Tng R = 200 
GeV. We want to stress that such a cancellation is not necessarily accidental, but could arise 
as a consequence of a particular alignment of the R- violating interactions in flavour space. The 
simultaneous presence of different couplings Ayj entails, in our basis, some flavour violation 
in the up sector. We expect therefore new effects in D°-D° mixing and in the decay modes 
D° — ► e + e~, D + — > 7r + e + e~. These processes at present do not set constraints on X' 1 j l more 
stringent than the one considered above We also remark that the bound in eq. (^j) can 
also be further relaxed by analogous cancellations among various A' m couplings with different 
indices k. 

In any case, the cl interpretation of the HERA data seems to suggest that B(K + — > -n^vv) 
is very close to its experimental bound, a prediction which can be tested in the near future 
in the ongoing Brookhaven experiment ||34|| . Notice however that, while m dfl determines the 
effective interaction responsible for K + — > 7i + uu, the HERA process is sensitive to rric L . In the 
absence of a complete theory describing all supersymmetric particle masses, m d - and are 
not necessarily related, and this prevents us from a definite prediction for B(K + — > 7i + uu). 

Notice that one mass relation can be obtained in the case of the cl interpretation of the 
HERA data, with the help of weak-SU(2) symmetry. Indeed, assuming no significant left-right 
squark mixing, we can predict m~ Sh = ^vn 2 ^ — cos2/3M^ ~ 200 to 220 GeV. The squark §l 
cannot be produced at HERA from valence parton processes, but could be observed at the 
Tevatron. 

The last possibility of a valence production mechanism is e + d — > ti via A' 131 . Apart from the 
HI and OPAL limits mentioned earlier, this coupling constant is constrained by experiments 
on parity violation in atomic physics, which imply |35 



A' 131 | < 4 x 10" 



m 



("20(1 GeV 



t L 



(10) 
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This certainly allows a sufficient production rate, even if B is significantly less than 1. In the 
absence of left-right stop mixing, one expects m- hh = \Jm\ L — cos2/3Mfy — mf + m| ~ 100 to 
130 GeV, in which case there would be an excessive contribution to the electroweak p parameter: 
Ap ~ 2 to 4 x 10~ 3 . At present the experimental value of Ap from LEP/SLD data plus mw/mz 
measurements is Ap exp = (4.7 ± 1.3) x 10~ 3 PS| . The Standard Model value for m t = 175 GeV 
is ApsM = 5.7 x 10~ 3 and 4.9 x 10~ 3 for m H = 100 and 300 GeV, respectively. Thus, there 
is little space for a new positive contribution to Ap. This suggests the necessity of significant 
left-right stop mixing, which is not unnatural and could also accommodate the lightness of the 
stop with respect to the other squarks. In this case, b\ could be heavier and Ap reduced. This 
also entails that the value of A' 131 inferred naively from the HERA data is smaller than the 
actual value by a factor of cos0(, where <pt is the stop mixing angle. 

Most sea production processes are excluded by a combination of different experimental 
constraints. Some of these have been discussed above; others come from contributions to the 
electron neutrino mass |57f 



a; 3 3I<5xio- 3 (-^-V , |A; 22 i<ixW-^r , en) 



,200 GeV/ ' 1 1221 V200 GeV 

which exclude sea-quark production mechanisms involving only third- (or only second-) gener- 



ation particles. In eq. (11) we have assumed a common supersymmetry-breaking mass wig in 
the 2x2 mass matrix for the d^-d^ system, and taken m s (mq) = 100 MeV for the running 
strange quark mass. Finally, ul production off sea quarks of the second or third generation is 



constrained by limits on charged- current universality, which impose 



|A - l<6Xir2 (2^fev) *»* = !, 2, 3. < 12 > 

Also, if the observed anomaly is due to the sea process e + u — *> dk R , then an effect more than 
50 times larger should have shown up in e~u — > dk R , while no anomaly has been observed in 
about 1 pb" 1 of data collected in e~p collisions PS 



The only remaining possibility for production on sea partons is e + s — > ti via the A' 132 
interaction, which is constrained, but not quite excluded, by the OPAL analysis ||11J| . Limits 
on anomalous top quark decay modes also set weak bounds on this coupling but these 
disappear as soon as the selectron mass is not much smaller than m t . This interaction can 
also give new contributions to the b — > decay rate, but these effects can be suppressed by 
an approximate alignment in the down sector. Therefore A' 132 could be as large as 0.3 and the 
process e + s — > ti be at the origin of the HERA signal. In this case, as we will discuss in the 
next section, the scattering process e + e~ — > ss has an anomalous contribution due to t-channel 
stop exchange which can be easily identified at future LEP 2 runs. 

4.2 Decay patterns 

Next we address the issue of the squark decay modes. In the case of cl, the most important 
possible decay modes are the i?-conserving cl — > cx® (i = 1, ••, 4) and Cl — > s\j (j = 1, 2), and 
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the R- violating c L — > de + , where Xi , Xj denote neutralinos and charginos, respectively. If im- 
parity violating couplings other than A' 121 were present, further decay modes could be allowed, 
although this possibility is severely constrained by limits on flavour and lepton conservation. 
The decay rate for the -R-parity violating mode is |39|, |40], |5T 



T{c L ^e + d) = ^-{\' 121 ) 2 m~ CL (13) 

1D7T 

and the coupling A' 121 is fixed by our production assumption. It has often been found that 
the .R-conserving modes dominate, but this is not necessarily the case. They could be either 
suppressed by phase space or, in the cl — > c\\ case, by (partial) cancellations in the neutralino 
couplings. The sxt decay mode can only be suppressed by phase space, so we assume that 
m + > 200 GeV. Neglecting cl r mixing, the decay rate for the neutralino mode is given by 

Xj 

0] 

a 2 ( m 2 o 

m c X °) = ^(A 2 + b 2 ) m , L . (14) 



32tt v \ ml 



where 



m N-a 1 

A * = U C B * = N i*+n t**0 w N a . (15) 

M w sm p 3 

In eq. fllBT) the iVy are the elements of the unitary matrix that diagonalises the neutralino 
mass matrix in the SU(2) - U(l) gaugino basis [P77| , and tan/? is the ratio of Higgs vacuum 



expectation values. The quark-mass-suppressed A 2 term in eq. (1141) may be neglected, and we 
notice that the B 2 term in eq. (fH]) is reduced either if the lightest neutralino is an approximate 
higgsino (Nu ~ N\2 ~ 0) or if there is a cancellation 

N 12 -tan^wr-^u • (16) 

Remarkably enough, we find that the cancellation ([16]) does occur in an acceptable domain of 
supersymmetric parameter space, given analytically by 

4sin 2 # iy 

m X? = o a ■ 2q M 2 I 7 

1 3 — 4 sm v w 



H = sin 2f3X ± ^/sin 2 2/3X 2 + m x o (m x? + 2X) (18) 

2(1 - sin 2 M(3 - 4 sin 2 M| 

(3-8sin 2 ^) M 2 ' 1 j 

Here M 2 is the SU(2) gaugino mass, while the U(l) gaugino mass is determined by the unifica- 
tion relation M\ = (5/3) tan 2 9wM 2 . Notice that such a cancellation is possible for c L — > ex®, 
but impossible, for instance, in the analogous d,R decay, whose rate is still given by eq. (^4|), 
with Bi = — tan 9w Nil. 

The results of numerical studies, see Fig. (|a, explicitly show the three regions where the 
cl -R-parity violating decay modes become important. The first region occurs for M 2 and /i 
large enough to suppress kinematically all two-body -R-parity conserving modes (m x o,m x ± > 
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Figure 4: Contours of B{ed) for the R-violating decay of cl, in the \i — M2 plane. 
Here A' has been fixed to 0.04 and tan/3 = 1 (a) and 5 (b), respectively. The LEP 2 
bound of 85 GeV for the chargino has also been implemented. 
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Figure 5: Contours of B(ed) for the R-violating decay oftL, in the 11 — M2 "plane. Here 
X' has been fixed to 0.0 Jj. and tan/3 =1 (a) and 5 (b), respectively. We have assumed 
a vanishing stop left-right mixing. The LEP 2 bound of 85 GeV for the chargino has 
also been implemented. 
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rric L ). In this case, B = 1 and di should lie at the edge of the parameter region excluded 
by DO. The second region is the thin slice of parameter space where x° is an approximate 
higgsino (fi << M 2 ). In this case, the couplings of cl to the light chargino and neutralinos are 
suppressed, and the .R-parity violating mode can compete with the i?-parity conserving ones. 
In the third region, the decay mode cl — > cx° is suppressed by the approximate cancellation 
of eq. (16). This cancellation is especially marked for tan/? ~ 1, where B > 0.5 over a large 
domain of /i. The extent of the cancellation region is reduced as tan/3 is increased, as can be 
seen from Fig. |]b, which is for tan/3 = 5, since the region where eq. ( |16|) is approximately 
satisfied becomes narrower as tan/? increases. We conclude from Fig. |] that the detection of 
e + q final states by HI and ZEUS data should not be a surprise. 

A small value of B could in principle be accommodated by increasing the magnitude of the 
A' 121 coupling by the corresponding factor of l/y/B, though the scope for this is severely limited 
by the bounds described above. On the other hand, if B ~ 1, the squark should be at the verge 
of being discovered at the Tevatron or possibly ruled out by CDF and DO data, as discussed in 
the previous section. 

In the case of fr,, it is interesting to notice that the neutralino decay mode II — > txl is 
kinematically closed in a natural way. In order to obtain a large value of B, it is sufficient to 
require that all charginos are heavy enough to forbid the decay £l — > bxf (see figs. [|a and 
||b). In view of the Tevatron limits discussed in the previous section, that may pose a problem 
if B(eq) is very close to 1, we have analyzed other possible decay modes. Under the gaugino 
unification assumption with m x + > m^ L , the three-body decay II — > bW + x^ (mediated by a 

virtual xf^t or bi) has a negligible rate. Even if a slepton were much lighter than the stop, 

the decay modes II — > 6£ + z>, bu£ + could not compete with the -R-parity violating decay, for 

^131 = 4 x 10~ 2 . It is usually assumed that the flavour-violating decay tt — > cx\ is rather 

suppressed in supersymmetric models. However, in theories with .R-parity violation, the whole 

issue of flavour conservation is undermined, and we cannot exclude new unexpected effects, 

which could lead to large rates for ti — > cx® and values of B considerably smaller than 1. 

If the stop is produced by the sea-parton collision e + s — > II, then the i?-parity violating 

decay is fast enough to compete with the chargino mode. The .R-violating branching ratio 

B depends only on the chargino masses m + and their gaugino compositions \Vji\. We find 

Xj ... 

B = 0.5 either for a pure gaugino-like chargino (|Vii| = 1) of 150 GeV, or for a mixed chargino 

(|Vii| = of 120 GeV. Therefore, in this case, it is easier to escape the Tevatron limits, 

although a small value of B requires a large value of A' 132 and a large effect in e + e" — ► ss at 

LEP 2. 

It is natural to ask whether the i?-violating scenario for the HERA events discussed above 
is compatible with the i?-violating scenario proposed elsewhere to interpret the four-jet 
excess seen by ALEPH [^3| at LEP 2 0. The suggestion contained in ref. |22|] was that ALEPH 
had observed the i?-conserving production process e + e~ — > eifiR-, mediated by an approximate 
U(l) gaugino with mass M\ between 100 and 120 GeV. The subsequent _R-violating decays of 
each e^/j into qq produce a pair of hadronic jets via an interaction of the A' type. This could 
have the 121 flavour structure advocated above for cl production at HERA, in which case we 

12 We are aware that the ALEPH signal [p3| has not been confirmed by the other LEP collaborations p4j , 
but reserve judgement on the final fate of the four-jet excess. 
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Figure 6: Crosss section contours for e~L&R ^nd vv "production at LEP 2 (yfs = 161 
GeV and tan/3 = 1). The calculations were performed using the programs SUSYXS 
documented in ref. /f^/, and include the effect of ISR. The regions corresponding to 
lightest chargino and neutralino masses heavier than 85 GeV lie above the dashed and 
dotted lines, respectively. 



predict the presence of a charm quark and an antiquark in the two dijets of the ALEPH final 
states. A strong constraint on this scenario comes from an adequate suppression of sneutrino 
pair production, which requires tan/5 ~ 1 and M<i larger than what predicted by gaugino 
unification. It is interesting to notice that such a choice of parameters can also lead to the 
approximate cancellation described in eq. (0), and therefore to a significant value for B. We 
display in Fig. |6| contours of the cross sections for e^eR and vv production at LEP 2 for y/s = 161 
GeV and tan/? = 1. We do not show the corresponding figure for larger values of tan/? since, 
for rrie L = 58 GeV, values of tan/3 > 1.23 are excluded by the requirement > M z /2. In 
Fig. |6] we have assumed gaugino mass unification, in order to allow the comparison with the 
results shown in Fig. |4j. However, as mentioned above, a better agreement with the ALEPH 
data is actually obtained when Mi is larger than (3/5) tan -2 9wM\. We conclude that the 
R- violating interpretation of the HERA data in terms of cl production advocated above is not 
incompatible with that proposed previously |22| for the ALEPH four-jet events. 
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5 Tests to Discriminate between Models 



In this final section we review some key experimental tests that may help to distinguish between 
different novel physics interpretations of the HERA large-Q 2 events. 

The first comment follows from Fig. [l]: the Q 2 distribution expected from effective contact 
interactions and resonance interpretations are different and could be distinguished clearly with 
a modest increase in statistics. The present data seem to favour a resonance interpretation, 
but it would be premature to draw firm conclusions at this stage. As for the x distributions, 
more statistics are again required to establish consistency with a resonant peak smeared out 
by ISR (see the Appendix), gluon radiation, hadronization and detector effects. 

"Charged current" events due to vq decays would be expected in some leptoquark and 
i?-violating squark scenarios at rates similar to those for the "neutral current" events seen. 
However, this is not the case, in particular, for the valence production e + d — > c^/t scenario 
favoured above. The HI collaboration has reported [|TJ four "charged current" events in a 
kinematic region where less than two are expected according to the Standard Model. The 
scenario we favour would be excluded if this signal built up into a significant signal with the 
advent of more data. 

The recorded luminosity in e~p collisions is rather limited, although it has already provided 
some constraints on scenarios in which a leptoquark (i?-violating squark) is produced via e + 
collisions with a sea quark, as commented in the previous section. A much higher e~p integrated 
luminosity should become available in the future. The cross section curves shown in Fig. |2] are 
also applicable to these collisions. It is a key prediction of our preferred ci/t model that there 
should be no large cross section for the production of a resonance peak in e~p collisions. 



One of the options for future HERA running is for e + D collisions El. This is potentially 



interesting, since the HERA signal is made from e + d collisions according to our favoured inter- 
pretation, the neutron contains twice as many d quarks as the proton, and this ratio is further 
enhanced at large x. The luminosity for e + d collisions in scattering off a neutron is the same as 
that for e + u collisions in a proton, and can be read off from Fig. 0. Unfortunately, the effective 
Eqm in e + n scattering will be less than the 300 GeV currently attained with protons. 

We recall that, for the reasons discussed above, the cl — > c(x° — > qq£, v) decay chain may 
have a branching ratio comparable to the e + q final state that we hypothesize to have been 
observed, and these should be observable at HERA. We note that £ ± final states are equally 
likely in \ decays, and that dominance by t = /i, r cannot be excluded. These final states would 
all have very clear signatures: charged leptons which may well have different flavour and/or 
charge from the incoming e + , or missing energy carried away by a neutrino, each accompanied 
by three hadronic jets, at least one of which should contain a charmed particle. 

In Figure we show the shapes of the x e Q and Q 2 distributions for wrong-sign leptons 
due to i?-conserving decays of squarks at HERA. We have not implemented any selection cut 
in preparing this figure, with the exception of a Q 2 > 5000 GeV requirement for the events 
in the x e distribution. Small Standard Model backgrounds to these final states are expected 
to come from semileptonic decays of heavy quarks (charm and bottom). The separation of 
these backgrounds depends significantly on the lepton isolation requirements, which are related 
to the detector characteristics and will not be studied here. Since the distributions depend 
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is the Bjorken x variable extracted using the electron method (see refs. jj], ||) 
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Figure 7: Kinematical distributions of wrong-sign electrons in the R-conserving de- 
cays of squarks, followed by R-violating neutralino decay. Different combinations of 
neutralino and slepton masses are shown. 



significantly on the masses of the states produced in the decay chain, we consider two values of 
the slepton mass (53 and 100 GeV) and three values of the neutralino mass (100, 150 and 180 
GeV). As expected, the figures show that the most interesting signals arise in the case of the 
largest neutralino masses and smaller slepton masses. 

Our analysis has pointed up the urgency of a joint analysis of the CDF and DO dielectron 
data, to see whether the existence of a 200 GeV leptoquark (or i?-violating squark) can be 
probed with the available Tevatron data, and, in the absence of a signal, down to what B(eq) 
it can be excluded. If the existence of such a leptoquark (squark) is still a live issue at the 
time of the next Tevatron run, we expect that the data taken then should be able to probe its 
existence down to values of B(eq) that are below those of interest to the present HERA data. 

The cascade decays of the squark via an .R-conserving cl —>■ X° c interaction could also have 
distinctive signatures at the Tevatron. We show in Fig. |8] the invariant mass distribution of 
lepton pairs, which is independent of the relative charge of the leptons. The final states consist 
of two leptons (with equal probability of having same or opposite charge), a large number 
of jets (there are 6 energetic quarks in the final state), and no missing transverse energy 

14 If there are also x ~~ * y v decays, more signatures would appear, such as e + jets + missing transverse 
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Figure 8: Invariant mass distributions of electron pairs from the R-conserving decays 
of squark pairs produced at the Tevatron, followed by R-violating neutralino decays. 
Different combinations of neutralino masses are shown. 



Possible backgrounds come from the production and decay of top quark pairs, when the net 
missing transverse energy carried away by neutrinos is small and additional jet are produced 
by radiative processes. In the case of like-charge leptons, the tt background requires one of the 
leptons to come from the semileptonic decay of one of the b quarks in the final state. Once again, 
the precise size of the background will strongly depend on the lepton isolation requirements, 
and will not be estimated here. 

In the case of stop production and decay via a virtual chargino, the signatures become 
particularly interesting, because of the presence of b jets in the final state. Transverse missing 
energy would also arise from the chargino decay to neutrino and slepton. In the case of cl decays 
to cx° there will be charm jets, which will still give rise to secondary vertices. Only higher 
statistics will however allow this signal to be isolated and distinguished from the secondary 
vertex distribution of b decays. 

It is natural to ask whether a squark weighing 200 GeV might have some observable indirect 
effects at LEP 2, even though it could not be produced directly. An R- violating interaction 
A' gives a contribution to the cross section of a generic quark- ant iquark final state // via the 

energy. 
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exchange of a squark q L or q R , which is parametrised as 



o = cr SM + 



3X' 4 h , 3\ ,2 a em l 



64ns 4s 



e f sjs - Mf ) 

e e e f + a ^L,R is _ M2)2 + T2M 



(20) 



where 



h = i±^-2x, l„fi±^) (21) 



1 ~t~ y 
1 

2 



I 2 = ± -*,- + *§ In (i±^£) (22) 



with a{ H = (T 3 ^ — e-^ sin 2 6 l H /)/(sin ^ cos^v^) and = — . In the case of dd production by 
cl exchange, we have a d R = — e^tan^, whereas in the case of cc production by d R exchange 
we have a c L = (1/2 — 2/3sin 2 6'vi/)/(sin6 l iycos6 l w). The contribution in eq.(pO|) proportional to 
Ii arises from the diagram with the R-parity violating vertices, while those proportional to I 2 
are interference terms with the Standard Model s-channel 7 and Z exchange respectively. For 
\fs = 192 GeV, the correction to the Standard Model cross section is ~ 0.02 pb, which we 
suspect that is rather small to be observed. However, there could be an observable signal in 
e + e~ — > ss due to t exchange in the sea production scenario discussed above, which is already 



on the verge of exclusion by OPAL [ I I 



Since the lightest neutralino x° decays rapidly via i?-violating interactions, the reaction 
e + e~ — > x°X° should be observable at LEP 2 for m x o below the kinematic limit, currently 
about 85 GeV. We do not discuss here the production cross section, which depends, e.g., on the 
selectron masses assumed. We have plotted in Fig. |6| the m x o = 85 GeV contours. Comparison 
with Figs. § and || indicates that our favoured i?-violating HERA scenarios are not strongly 
constrained by present LEP 2 data, though future data at y/s = 200 GeV might make some 
inroads on the parameter space. 

The .R-violating scenario squark mentioned in the previous section may have observable 
consequences for other experiments that are not a priori related. One example is K — > -nvv 
decay pi ]. We have seen that this imposes one of the most stringent constraints on the A' 121 



coupling that we invoke. A corollary is that there may be an interesting contribution to this 
decay from beyond the Standard Model, waiting to be discovered just below the present level 
of experimental sensitivity. The magnitude of any such signal depends on the pattern of flavour 
mixing among i?-violating couplings. In some variations, there may also be contributions to 
nuclear (3(3 decay lurking just below the present experimental sensitivity p0| . 

These examples point to a general theoretical issue raised by the possibility of R violation. 
General /^-violating couplings do not respect the classic conditions for natural conservation of 
flavour in neutral interactions ||44|| . Perhaps these are in any case optional, and one should 



be content with models which fall numerically below the experimental upper limits on flavour- 
changing neutral interactions. On the other hand, natural respect for them played an important 
historical role in motivating the Standard Model. Therefore, it is desirable to clarify whether 
there are any interesting and plausible conditions under which these constraints are naturally 
respected by R- violating interactions [^, [RJ. 
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This example shows that interesting and relevant theoretical, as well as experimental, issues 
are into new light by the observation of large-Q 2 events at HERA. It may well be that these 
turn out to be a malign statistical fluctuation, rather than a harbinger of new physics. However, 
we have shown in this paper that complementary experiments may soon be able to cast light 
on possible interpretations in terms of physics beyond the Standard Model. In the mean time, 
the HERA large-Q 2 events have caused us to look anew at supersymmetry with R violation, 
in particular, and provided us with new reason to question the conventional i?-conserving 
paradigm for supersymmetric phenomenology. 
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Note Added 



After the completion of this paper we received the articles in ref. [(47]> fUL which make the 



interesting observation that experiments on atomic parity violation pose stringent limits on the 
strength of any A e V q = RR — LL + RL — LR eeqq contact term. The present experimental value 
on the coefficient C\ q of the A e V q term in the effective lagrangian, as given on pages 87-92 of 
ref. ||25|| , implies the following 95% CL limit on the corresponding deviation from the Standard 



Model AC 



1<2 



Y^ij(V27rr)ij5i) / (GF-hfj) (where 5i = +1, — 1 for i = R or i = 
- 0.099 < AC lu < 0.051 - 0.050 < AC ld < 0.084 



Assuming a single non-vanishing operator at a time, this translates into 



AH A^ > 2.0 TeV, 



A LL> A Li?> A i?L> 1 ^RR >z -< iev 5 

A+f, Att A^?, A D t>2.8TeV, 



RR 



Alt A7i, A+£, A+£ > 2.1 TeV. 



L, respectively): 
(23) 

(24) 

(25) 
(26) 
(27) 



Comparing with Fig. |T] we see that these limits are quite constraining on the individual terms. 
But we also see from Fig. H that we could, for example, take Aj^ — Aj^ [0, which is parity 



conserving, and add the corresponding contributions that are very similar in shape. In this way 
the limits are evaded and the fit is as good as for the separate contributions. 



Appendix: Some Attempts to Understand the Possible 
Effects of Initial-State Radiation 

In an attempt to gain more insight into the apparent spread in the invariant masses of the ZEUS 
events, and the fact that their masses appear at first sight to be somewhat higher than those 
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of the HI events, we have looked into the possible impact of initial-state radiation (ISR) In 
the presence of ISR, the relation between the value M e of the reconstructed resonance mass, as 
determined by the electron method, is related to the true value M by 

Mj = M 2 |^-p (28) 

where z is the fraction of the electron's longitudinal momentum lost to ISR, and y e is the value 
of the conventional deep-inelastic y variable estimated using the electron method. In the case 
of small z, eq. reduces to: 



Ml = M 



l + z 



' 22/ e -r 

. Ve , 



(29) 



which corresponds to a negative shift in mass for y e < 1/2, and to a positive shift for y e > 1/2. 

The analogous relation between M 2a , the mass determined by the double-angle method, 
and the true value M is 

M L = M2 ^hf ■ (30) 

In the presence of ISR, M 2a will therefore always be larger than the true value of M. In 
particular, M 2a will always be larger than M e . 

We recall that HI prefers to estimate M using M e |l|, whereas ZEUS favours M 2a @- ISR 
effects could therefore lead qualitatively to Mzeus being greater than Mhi, as observed, and 
resolution differences might explain their greater spread. On the other hand, the experimental 
cuts allow for a fraction of electron energy lost to ISR up to ~ 10%, so it is not clear whether 
its effect could be important quantitatively. 



One can use eqs. ([280 an d (30) to extract a relation between the masses reconstructed with 
the two techniques and the true mass, under the hypothetical assumption that the measured 
differences are due to ISR and not to resolution effects. The following relations hold: 

M 2 

z = y e (l-p) , P=j7r (31) 

1V1 2a 

M 2 = M 2 a (l-y e + y e p) 2 . (32) 

We applied these relations to the five ZEUS candidate events, using the values of y extracted 
using the double-angle techique. Four out of the five events have x e < x 2a , and are therefore 
compatible with the ISR interpretation. The values of z and of the true mass which we calculate 
using the above relations are given in Table [| 

Allowing either the highest- or lowest-mass event to be background, we find good consistency 
with a single mass value around 220 GeV, though this still looks higher than that quoted by 

hi m 

15 It is clear that this exercise is best carried out by the experimental collaborations themselves: our intention 
is only to form an approximate impression of how large the effects of ISR might be. 

16 We have checked that, when applied to the HI data, the above ISR estimation procedure makes no significant 
difference to their preferred mass value M ~ 200 GeV, which has an energy scale error of about 5 GeV. 
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Table 3: Estimated effect of ISR on the apparent masses of the five candidate ZEUS 11/ 
events extracted using the electron (M e ) and the double-angle (M% a ) techniques, where 
z is the energy fraction lost to ISR as estimated from eq. $3\ ) and M is the ISR- 
corrected mass defined in eq. ffiy- The first event gives no indication that z > 0, and 
we have retained the M.2 a estimate. 



Needless to say, only an accurate analysis by ZEUS, properly accounting for the effects of 
experimental resolution and their correlations between the two techniques, will provide an ac- 
curate estimate of the ISR- induced corrections. As remarked in ref. 0, the differences between 
the values of M e and M 2a observed in the five candidate events are also not inconsistent with 
the reported measurement uncertainties. 
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